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ABSTRACT 

The evolutionary history of many species exhibits episodes of habitat expansions and 
contractions, often caused by environmental changes, such as glacial cycles. These range 
changes affect the dynamics of biological evolution in multiple ways. Recent microbial 
experiments suggest that enhanced genetic drift at the frontier of a two-dimensional 
range expansion can cause genetic sectoring patterns with fractal domain boundaries. 
Here, we propose and analyze a simple model of asexual biological evolution at ex- 
panding frontiers to explain these neutral patterns and predict the effect of natural 
selection. Our model attributes the observed gradual decrease in the number of sectors 
at the leading edge to an unbiased random walk of sector boundaries. The long time 
sectoring pattern depends on the geometry of the frontier. Whereas planar fronts are 
ultimately dominated by only one sector, circular colonies permit the coexistence of 
multiple sectors, whose number is proportional, in the simplest case, to the square root 
of the radius of the initial habitat. Natural selection introduces a deterministic bias in 
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the wandering of domain boundaries that renders beneficial mutations more likely to 
escape genetic drift and become established in a sector. We find that the opening angle 
of those sectors and the rate at which they become established depend sensitively on 
the selective advantage of the mutants. Deleterious mutations, on the other hand, are 
not able to establish a sector permanently. They can, however, temporarily "surf" on 
the population front, and thereby reach unusually high frequencies. As a consequence, 
expanding frontiers are susceptible to deleterious mutations as revealed by the high 
fraction of mutants at mutation-selection balance. Numerically, we also determine the 
condition at which the wild type is lost in favor of deleterious mutants (genetic melt- 
down) at a growing front. Our prediction for this error threshold differs qualitatively 
from existing well-mixed theories, and sets tight constraints on sustainable mutation 
rates for populations that undergo frequent range expansions. 



1. Introduction 



Population expansions in space are common events i n the evolutionary his t ory of many species 
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from where they first evolved, invade into favorable habitats, or move in response to environmen- 
tal changes, such as the recent climate warming, glacial cycles, or gradients in nutrients, salinity, 
ambient temperature, etc., in the case of biofllms. Some species undergo range expansions rarely, 
because environments change slowly, others like epidemic pathogens do so frequently as part of 
their ecology. 

These range expansions cause strong differences between the genetic diversity of the ancestral 
and the newly colonized regions, because the gene pool for the new habitat is provided only by 
a small number of individuals, which happen to arrive in the unexplored territory first. The 
associated alteration of the gene pool depends on the specific demographic scenario, and encodes 
precious information about the migrational history of a spec ies. These "gene tic footprints" offer 
ways to infer, for i nstance, how humans moved out of Africa (jTempletonl I2OO2I ) or species respond 
to climate change (jHewittlboOol ). However, the underlying question, how to decip her the observed 
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The most widely appreciated consequence of a range expansion is a genetic bottleneck. Newly 
colonized regions are founded by a small subset of a larger ancestral population, typically with 
a greater genetic diversity. Because of this moving bottleneck at the expanding frontier, one 
expects a spatial gradient in the genetic diversity indicating the expansion direction of the ancestral 
population. The magnitude of this gradient depends sensitively on the population dynamics of 
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the pioneers at the frontier (e.g., Allee effects (1 Ailed 119311)). but only weakly on the maximum 
population density, also known as carrying capacity (IHallatschek and Nelson! 120081 ). 



These predicted gradients or "clines" in genetic diversity h ave indeed been picked up by large 
scale genomic surveys across populations (IHandley et al.1 120071 ). and provide valuable information 
about the demographic and ecological history of species. For instance, a freq uently observe d south- 
north gradient in genetic diversity ("southern richness to northern purity" (jHewit till 9961 )) on the 
northern hemisphere is thought to reflect the range expansions induced by the most recent glacial 
retreat. In the case of humans, the genetic diversity decreases essentially linearly wit h increas- 
ing geographic distance from Africa ([Rosenberg et al.l 120031 ; iRamachandran et al.ll2005l ). which is 
thought to be indicative of the human migration out of Africa. 

Until now, inference techniques have not made use of spatial correlations in the direction trans- 
verse to the gradient of genetic diversity. There is recent evidence from a microbial experiment, 
however, that these spatial correlations might actually be quite pronounced. Microbial systems have 
been establi shed over the last decad e as a valuable tool to probe fundamental aspects of evolution- 
ary biology (jElena and Lenskill2003l ). While microbial evolution experiments were fir st designed for 



well- mixed populations, they are now extended to spatially structured populations (IKeymer et al 



20061 ) . With these spatial sy stems, the genetic impa ct of range expansion on the genetic diversity 
has recently been measured (IHallatschek et al.ll2007l ). It was found that range expansion leads to a 
striking population differentiation along the frontier of the advancing front. As a consequence, sec- 
toring patterns emerge as a distinct footprint of past range expansions, see Fig. [TJ These patterns 
appeared, after two differently labeled, but otherwise identical, sub-populations were mixed and 
plated on a Petri dish. The initial liquid deposition on the Petri dishes took the form of circular or 
linear droplets, the later being inoculated off a sterile razor blade. As these initially well-mixed pop- 
ulations grew colonies, the mutant strains segregated at the wave front and gave rise to a sectoring 
pattern. The growth rate and cell mobility decrease markedly once the population wave passes by, 
leaving a frozen record of gene segregation in its wake. Qualitatively similar patterns were found in 
two very different microbial species, the bacterium E. coli (a mutant strain which lacks flagella) and 
brewer's yeast, Saccharomyces cerevisiae in its haploid form. The generi c nature of the sectoring 
mechanism sugg ests that sectors could be widespread in wild populations (lExcoffler and Ravi 1 2008: 
Biek et al. 2007 ). although direct evidence for those patterns is limited so far fjBiek et al. 2007 ). 
Detecting the trace of such spatial correlations in the genetic structure of a species could reveal 
details of ancient migration patterns. 

Although one motivation of this paper is a better understanding of populatio n genetics during 
range expansions, another is to explore the use of experiments such as those in Ref. (IHallatschek et al 
20071 ) as an assay to measure the effects of beneficial and deleterious mutations. Imagine mixing 
together a stable "background" or wild type strain of bacteria or yeast, labeled by a constitutively 
expressed green florescent protein with, say, a small population (e.g., 2 — 5% by cell number) of 
mutants, all labeled red. For simplicity, assume that all favorable mutants have an identical fitness 
advantage over the wild type, and that there is a similar identical fitness detriment for the dele- 



Fig. 1. — Competition experiments between neutral strains in growing bacterial (E. coli) and yeast 
colonies (S. cerevisiae). These colonies were grown from mixtures of CFP (red) and YFP (green) 
label ed cells, which were dep osited on the agar plates as linear (a,b) or circular (c,d) droplets (see 
Ref. flHallatschek et all 120071 ) for experimental details). Even though both strains of each species 
were otherwise genetically identical, the growing colonies exhibit a striking segregation of the two 
neutral markers (CFP and YFP) over time. The dynamics of segregation is restricted to the edges 
of the colony, while, except for a gradual thickening, the interior distribution of CFP and YFP is 
essentially frozen. Upon comparing the bacterial with yeast colonies (a,b versus c,d), it is apparent 
that the wandering of domain boundaries, and hence the genetic drift, is more pronounced for 
E. coli than for S. cerevisiae. In Figures c) and d), we have indicated quantities relevant to our 
quantitative analysis of the observed sectoring. The dashed box and ring in (c) and (d), respectively, 
indicate the initial size extent of the founding population. The size of domains at the expanding 
fronts are parametrized by the linear separation X(r) as a function of the effective "time" r in the 
linear inoculation, and by the sectoring angle <3>(r) as a function of radius r in a circular colony. 
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terious mutants. A possible outcome of a linear inoculation from liquid culture on a Petri dish is 
indicated schematically in Fig. As discussed in the Figure caption, this setup allows statistics 
to be gathered from numerous sectoring events in the same experiment. Beneficial mutations give 
rise to sectors with a characteristic opening angle <3>, although measurements could be obscured by 
genetic drift at the frontier. Deleterious mutations, which would die out rapidly in a completely 
deterministic scenario, can nevertheless take advantage of genetic drift to "surf" for a period of 
time. Very approximately, one can think of the march of the two population waves away from 
the frontier on the Petri dish as being like serial dilution experiments in liquid culture, with the 
distance from the homeland roughly proportional to the number of repeated "dilutions" . From this 
point of view, advancing population fronts are a low tech, massively parallel form of serial dilution, 
and could thus be used for evolution experiments. 

In this article, we develop a quantitative model of biological evolution at expanding frontiers, 
inspired by the microbiological experiments described above. Specifically, we focus on the sectoring 
dynamics after well-defined monochromatic domains have been established in the early stages of 
the experiment. Then, we argue that allele frequencies change due to the growth and shrinkage of 
sectors, which in turn reflects the competition of the deterministic force of natural selection and 
random number fluctuations (genetic drift). 

Even in cases where the temporal variation of sector sizes is purely stochastic (no selection), a 
gradual decrease in the number of sectors is expected as the population expands. This coarsening 
process occurs because domains frequently lose contact to the wave front and no longer participate 
in the colonization process Fig. [5l In Sec. 13.11 we treat this form of neutral evolution using a 
simple model of annihilating random walkers, which allows us to predict the gradual decrease in 
the number of sectors as the colony grows. Both linear and radial inoculations are considered. 
Although our model predicts that one allele (genetic variant) dominates for linear inoculations, the 
number of surviving sect ors remains finite at lo ng; times for radial expansions, in agreement with 



the experiments in Ref. (IHallatschek et al.l 120071 ). In Sec. 13. 2\ we consider biased random walks to 
study the spread and ultimate fate of beneficial mutations arising at the front of a population, such 
as the one documented in Fig. [3l We relate the asymptotic sector angle to fitness differences and 
determine the rate at which beneficial mutations become established based on their mutation rate. 
Finally, we determine the genetic load due to deleterious mutations that accumulate in the course 
of a range expansion. Our analysis shows that, due to enhanced genetic drift, selection is quite 
ineffective in purging deleterious mutations from the invasion front. Finally, using a combination of 
theoretical arguments and computer simulations, we find that a critical mutation rate exists along 
a linear frontier beyond which the front population would inevitably decline in fitness. Depending 
on the details of the range expansion, our prediction for this error threshold can be much lower 
than the well-known result for well-mixed "zero-dimensional" populations. 

The emergence of well-defined domain boundaries, which underlies our analysis, is not specific 
to range expansions in microbial systems, as we argue in Sec. HI which contains conclusions and 
discussion. Instead, sharp boundaries appear naturally due to the reduction of dimensionality (2 to 
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1) at the advancing front of a population spreading across a surface. Building on this hypothesis, 
we discuss to what extent our results generalize beyond microbial populations. 



2. Materials and methods 

Our null model for the sectoring dynamics rests on the assumption that the reproductive 
success of an individual is independent of its color. In other words, individuals with different 
alleles (genetic variants) have the same fitness. In a strictly neutral population mixture, domains 
should not grow on average, except for a line ar increase in the radial direction due to the inflating 



frontier in the case of circular inoculations (jHallatschek et al.l 120071 ) . Even in the neutral case 



however, sectors will exhibit variations in their sizes, because of inevitable chance effects during the 
reproduction process (genetic drift). The variations in sector sizes are manifest in the erratic path of 
the sector boundaries, see Figs.[]Ji,b. Apparently, these effects are much more pronounced for E. coli 
populations than for S. cerevisiae. Due to these fluctuations, neighboring domain boundaries can 
collide as the colony grows larger. As a result, the enclosed domain loses contact to the colonization 
frontier and is henceforth trapped in the bulk of the colony without further participation in the 
colonization process. Thus random fluctuations due to genetic drift are responsible for the continual 
reduction in sector number. 

Our phenomenological model for these neutral dynamics is illustrated for a linear inoculation 
in Fig. [H The meandering ends of domain boundaries are represented by random walkers that 
populate the expanding edge of the colony. The random walkers come in pairs, move along the 
advancing frontier and annihilate when they meet. Their trajectories in space and time describe 
the meandering and coalescing of domain boundaries visible in the microbial experiments of Fig. [TJ 

This model of annihilating random walkers with dynamics embodied in domain walls instead 
of the cells themselves serves as an effective description of the neutral evolutionary dynamics on 
long time and length scales (see discussion). The mathematical formulation of the model leads to 
generalized diffusion equations, which could be solved analytically. 

The evolutionary dynamics changes dramatically when mutations arise close to the range 
margins of an expanding population that have non-negligible fitness effects. The funnel-like sector 
shape on the left side of Fig. 02 for example, is the result of a mutation that increased the rate of 
expansion in this particular green population compared to either its unmutated green relatives or 
the red population. A series of similar observations shows that beneficial mutation generically give 
rise to sectors with unusually large sector angles. This phenotype implies a bias in the diffusion 
process of the sector boundaries. To account effects of selection, we thus formulated a biased 
diffusion model for the domain boundary motion. This model was solved analytically in the limit 
of small mutation rates, and studied by means of spatially explicit simulations for larger occurrence 
rates of deleterious mutations. 
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Fig. 2. — Schematic illustration of the result of a linear inoculation along a Petri dish of a small 
percentage (say, 2 — 5%) of mutant cells with both beneficial and deleterious mutants (labeled red, R 
and R' respectively) with a majority of wild type (labeled green, G). (The width of the "homeland" 
between the two dashed lines is not relevant for this discussion, but is determined by the amount of 
fluid on the razor blade and the capillary length of the liquid culture). As discussed in the text, the 
opening angle $ (shown for a beneficial mutation at the bottom) is given by the square root of the 
fitness advantage. The potentially obscuring effects of genetic drift in the linear boundaries of the 
beneficial mutant strains are not shown. Two deleterious mutants that have "surfed" successfully 
for a time at the population front are shown as well. These surfing events owe their very existence 
to genetic drift. 




Fig. 3. — Beneficial mutations give rise to sectors with unusually large opening angles, a) This 
colony of yeast (S. cerevisiae) was grown from a 5:1 mixture of CFP (red) and RFP (green) labeled 
cells. The large funnel-like green sector, which arose spontaneously, outgrows both wild type strains, 
b) A linear inoculation of mixture of a CFP (red) strain that has a beneficial mutation compared 
to an otherwise neutral RFP (green) strain. The resulting sectors of CFP mutants have similar 
shape and are well separated due to the small ratio (1:40) of mutant to wild type. 
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3. Results 

3.1. Neutral evolution 

We first present scaling analysis for the emerging neutral coarsening dynamics, at both planar 
and curved fronts. The results of a more precise analytical calculation of the sectoring dynamics 
follows thereafter. 



3.1.1. Scaling analysis 

Consider the dynamics of the distance X(r) between the tips of two neighboring domain walls 
within a linear frontier, such as in Fig. We assume that the domain boundary separation 
X(r) is a continuous random variable HI that fluctuates as if each domain boundary carried out an 
independent random walk with diffusion constant Dx- This assumption implies that if the average 
front position advances from ro to r (see Fig. [Tfc) by a length increment Ar = r — ro, then the 
associated increment AX = X(r) — X(tq) in the domain size has zero mean and a variance a\ 
that grows linearly with distance Ar, 

(AX) = , (AX 2 ) = 4D x Ar = <j\ • (1) 

Here, angular brackets denote an average over many realizations and the diffusion constant Dx 
describes the wandering of a single wall and has units of length 2 /length. The extra factor of 2 
arises in Eq. (pQ) because we look at the difference coordinate between two independent random 
walks [E Note that Eq. (pQ) also implies a random walk in real time. Indeed, if the colony expands 
at a constant velocity v, we have Ar = vAt and may write Eq. (pQ) equivalently as 

<7x = 4:DxvAt = ADxAt . (constant expansion speed v) (2) 

Thus, the random variable X carries out a random walk in real time t with diffusion constant 
Dx = DxV- In the following, we will employ r as the effective time-like variable rather than the 
real time t. Although the expansion velocity v will vary from organi sm to organis m, and may even 



be time dependent during the early stages of a radial expansion (jMurrayl 12004 ) , v(t) drops out 
when the problem is formulated as in Eq. (pQ). 

The random walk assumption, Eq. (ED) , forms the basic work ing hypothesis of our analysis 



of neutral evolution. As discussed in Ref. (jHallatschek et al.1 120071 ). this assumption can be vio 



lated when the interface of the expanding population is rough, which can drive a super-diffusive 



1 Throughout this paper, we denote random variables by capital letters and their values for specific realizations by 
corresponding lower case letters. 

2 If the wandering of neighboring domain boundaries is correlated, e.g., due to interactions between the walls, we 
may consider Eq. JT|) simply as a definition of the phenomenological parameter Dx(r). 
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wandering. This super-diffusion can also be analyzed on the scaling level, and we shall indicate 
the implied changes in Sec. [H The diffusive scaling used here nevertheless captures the essential 
behavior of many models and has the advantage of being exactly solvable. Note that our neglect 
of the roughness of the interface allows us to characterize the position of the population front by a 
single time-like variable r for both linear and circular inoculations (see Fig. [1]). 

A coalescence event occurs when two domain boundaries meet, and the enclosed domain is 
displaced from the wave front in favor of the neighboring domains. Equivalently, we can say 
that the tips of domain boundaries annihilate when they meet. The typical distance between 
domain walls that evade annihilation at distance r should be comparable to the standard deviation 
dj(r). Hence, we show that the average number N(r\ro) of surviving lineages at effective time r 
declines with increasing Ar as the inverse of ax- Using a pre-factor derived below in the diffusion 
approximation, the number of sectors then is given by 

N(r\r ) = — = J == . , (linear inoculation) (3) 
V27R7X 2y/27rD x Ar 

for a neutral 50 : 50 mixture of two differently labeled but neutral sub-populations. In Eq. ([3|), L 
is the total length of the population front. Note that microscopic length scales, such as the cell 
size (or, more generally, an initial domain size) do not appear in this asymptotic formula. The 
annihilation process eventually leads to the prevalence of only one domain, or allele, after the front 
advances a distance Ar such that the distance between domain boundaries has become of order L. 
Although Eq. ([3j) is only approximate in this regime (see below), it correctly suggests Ar ~ L 2 /Dx 
as the order of magnitude of the fixation time for a linear inoculation. 

In the case of an advancing curved population front, the ultimate fate of the gene pool is 
rather different. On top of the diffusion process, neighboring domain boundaries are subject to 
a deterministic expansion caused by the changing size of the perimeter of the population. Thus, 
diffusion competes with an antagonistic drift caused by the inflation of the perimeter, which inflates 
the distance between neighboring domain walls. Although the deterministic drift term will dominate 
on long times, diffusion dominates at earlier times. 

Let us consider a circular colony as the simplest example of an advancing curved front. In 
this case, the perimeter grows linearly with the radius of the colony. This radial expansion tends 
to increase arc-length distances between domain boundaries. The ends of two neighboring domain 
boundaries drift apart at a speed proportional to their distance. Thus, two neighboring random 
walkers are subject to both deterministic drift and diffusion in the arc-length parametrization of 
their separation. Alternatively, we can describe the distance between two neighboring random 
walkers by their angular distance 

A$ EE — , (4) 

r 

as indicated in Fig. [Hi. This change of variables simplifies the problem to pure diffusion because 
sector angles remain unchanged on average during the growth of the colony. The angular diffusion 
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constant D$ becomes, however, a decaying function of the radius r of the colony, 

(A$ 2 ) (AX 2 ) jr 2 D x 

D*(r) = lim ^—^ = lim V J 1 = , (5) 
v 7 Ar^o 4Ar Ar^o 4Ar r 2 

where is the diffusion constant that describes wall wandering for a linear inoculation, Eq. (pQ). 
When the colony grows from an initial radius ro to a final radius r, the mean square angular 
displacement changes by 

4 = <A$ 2 > = 4 fdr' D<s>{r f ) = 4D X (r^ 1 - r' 1 ) . (6) 

Notice that the variance in angular distance depends on increments in inverse radii, as opposed to 
the linear dependence Eq. (pQ) for arc-length distances. An immediate consequence is that, when 
the colony grows much larger than the initial radius, r/ro ^> 1, the expected mean square angular 
displacement stays finite, cr$ — > 2^/ Dx /r^. At long times, the sectoring pattern should therefore 
decompose into a finite number Af(oo|ro) given by 



7V(oo|ro) = — = • (circular inoculation) (7) 

The numerical pre-factor in Eq. Q again results from the quantitative analysis of the sectoring 
statistics in the diffusion approximation discussed below. 



3.1.2. Sector size distribution for linear and radial expansions 

According to our model, we can view successful sectors that do not get trapped behind the front 
as being bounded by random walks that evade any collision. The size distribution of sectors should 
therefore be determined by the positional distribution function of annihilating pairs of random 
walkers conditional on survival. Our goal is to determine this distribution quantitatively. 

As in the scaling discussion, we measure time by the spatial position r(t) of the frontier of the 
population. In the scaling discussion, we described the size of a sector by a distance X(r) in the 
linear inoculation, and by an angle $(r) in the circular case. In order to capture both scenarios 
simultaneously, let us introduce a generalized sector size variable Z(r), which is assumed to carry 
out a random walk with diffusion constant D{r), 

(AZ) = , (AZ 2 ) = D(r)Ar . (8) 

The specific scenarios of linear and radial inoculations can be recovered by identifying 



linear inoculation: Z = X , D(r) — Dx — const, 
circular inoculation: Z — $ , D(r) = D$(r) = Dx/r 2 



(9) 
(10) 
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The statistical properties of the sector size Z(r) are described by a diffusion equation. Let 
F(z, r\zo, ro) be the probability that Z(r) — z when the front is at r given that Z(ro) — zq at the 
earlier front position ro- This distribution function satisfies 



d r F(z, r\z ,r ) = 2D(r) d 2 z F(z, r\z , r 



(11) 



which is a direct consequence of Eq. (JSJ) and the continuous nature of the random variable Z(r) (jvan Kampen 

To account for the possibility of annihilation, we impose an absorbing boundary condition 
at z = 0, 

F(0,r|z ,r ) = 0. (12) 

The factor of two in front of the diffusion constant in Eq. (fTTI) arises because Z{r) is the distance 
between two random walkers. Also note that, in order to keep the analysis simple, we assume that 
space is unbounded in this section. Finite size effects can lead to the important effect of fixation 
as discussed in Sec. 13.1.31 

The absorbing boundary condition Eq. (fT2l) can be fulfilled exactly by writing the solution as 



F(z, r\z Q , r ) = G(z, r\z , r ) - G(-z, r\z , r ) 



(13) 



in terms of the solution G(z, r\z$, ro) of the diffusion equation (fTTI) without annihilation. Diffusion 
without annihilation, on the other hand, is well-known to be described by a Gaussian probability 
distribution, 



g(Az, Ar) = G(z + Az, r + Ar\z 0l r ) 



(14) 



V27TCT 2 

where a is the standard deviation of the random variable Z accumulated from frontier distance ro 
to distance r, 

2/ \-A[ r jir«i\ / 4D x (r-r ) linear 
a (r, r ) = A jjr D(r ) = | _ ^ ^ (15) 

This quantity was evaluated in the scaling analysis in Eqs. (pQ) and © for the linear and circular 
inoculation, respectively. Upon combining Eqs. (|T3l) and (TT4]) . we obtain 



F(z + Az, r + Ar\z 0l r ) 



g(Az, Ar) - g(-2z - Az, Ar) 



exp 



Az 2 
2a 2 



V2 



7TCT Z 



exp 



-2z (zo + Az) 



(16) 



This result can be used to predict the sector size distribution that emerges when a linear or circular 
colony grows from a finely divided mixture of differently labeled sub-populations. Suppose that 
the front of the colony advances from position ro to position r > ro. Then, every point along the 
edge of the population sends out domain boundaries in a tree-like web that gradually coalesces, as 
illustrated in Fig. El Due to this coarsening process, the number of sectors N(r\ro) at the front is 
a decreasing function o f time r, whose ensem ble average we seek to determine. Consider first an 
"infinite alleles model" fjHartl and Clarklll997l ) where each cell of the founder population is labeled 



a A h 

average I D 

growth | ?C~^ 

direction | ^ B 

Fig. 4. — An illustration of the mechanism by which the sectoring pattern of a microbial colony 
coarsens in time, (a) Four monochromatic domains are bounded by a moving frontier (black 
line), (b) As the colony grows further in the upward vertical direction, the domain boundaries 
follow wandering paths. By chance, the two domain boundaries on the left-hand side meet. As 
a consequence, the enclosed domain (red left-hand domain in (b)) loses contact to the population 
front and is, henceforth, trapped in the bulk of the colony. We model these dynamics by replacing 
the tips of the domain boundaries by random walkers (blue circles) that "live" on the growing one 
dimensional edge of the colony. Even though these random walkers annihilate when they meet 
(yellow star), they have a non-zero survival probability on the growing circumference of a circular 
colony. 




Fig. 5. — A sketch of the sectoring pattern emerging from an (hypothetical) initial population 
in which each individual is labeled differently. The curved colony grows from initial radius ro to 
radius r. Each surviving lineage emerges from a single founder cell denoted by a cross. 
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differently. In this case, each surviving sector at r originates at earlier effective time ro in a single 
individual, the most recent common ancestor (denoted by a cross in Fig. [5j). 

Consequently, a sector is generated by two domain boundaries that start from the same location 
of the common ancestor at initial time ro- Hence, the sector size distribution is simply given by 
the probability distribution function of random walkers that evade annihilation given that they 
start (almost) at the same place. This distribution can be obtained by normalizing the small zq 
expansion of F(z, r\zo, ro) in Eq. (fT6l) . We thus obtain the sector size distribution 

P(z,r\r ) = ^exp(-^j , z>0. (17) 

P(z,r|ro) is the probability that a randomly chosen sector sampled at time r has size z in the 
infinite alleles model. The product aP is plotted as a function of the dimensionless variable z/cr 
in Fig. [S] (dashed line). Note that Eq. ffTTD assumes that the population front is unbounded. For 
a circular inoculation, < z < 2tv, Eq. (fT71) is only valid in the limit a(r, ro) <C 1. The more 
complicated exact distribution for this bounded case, Eq. ([391) , is derived in Sec. 13.1.51 below. 

The mean of the distribution Eq. (fTTl) is given by 

\^(r,r ) (18) 



^2nDx{r — ro) (linear inoculation) 



'2ttDx (r 1 — r x ) (circular inoculation) , 

and represents the average size of a sector. Equivalently, it is the average distance between two 
non-colliding random walkers that initially start out from (almost) the same place. The numerical 
pre-factor may be compared with the expected size cj^/2/ti of a sector in the absence of annihilation, 
as can easily be shown from Eq. (fT4l) . Thus, the average separation of surviving annihilating random 
walks is a factor of tt/2 w 1.57 larger than of "phantom" random walks that can pass freely through 
each other. Intuitively, this factor represents an effective repulsion between two random walks that 
must avoid collision to survive. 

The average number 7V(r|ro) of sectors measured in many repeated experiments of a given 
kind (circular or linear), is given by the ratio of total size of the front L and average sector size, 

L [~2 L 

N(r\ro) = | — yj = Y (infinite alleles model) . (19) 

Here, L would be the total length of the population front for a linear inoculation. In the circular 
case, however, one would choose L = 2tv. The standard deviation then approaches a<$> = 2^/Dx/ro 
for r — > oc, so that 7V(oo|ro) = ^2nDx /tq. H 



3 Provided that sector interactions (neglected here) are short range, we expect that corrections to Eq. ([T9]) are of 
order y/N(r\ro). 
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These results, however, only hold in the infinite alleles model. When we relax the assumption 
that every founder individual has a different color, the average number of sectors will be less. For 
instance, if the init ial population is labele d by only two colors in equal proportions, as in the 
experiment in Ref. (jHallatschek et al.1 120071 ), two neighboring sectors of the infinite alleles model 
will have the same color with probability 1/2. Thus, there will be half as many sector boundaries 
as in the infinite alleles model. Accordingly, the average number of sectors in the two alleles model 
will be given by N(r\ro) = L(r)/ (2Z(r|ro)) 3 which gives the pre-factors used in Eqs. ([3j) and 
(|7|) of our scaling analysis. More generally, we may consider an initially well-mixed population, 
in which two r andomly chosen indiv iduals have a different color with probability H, known as 
heterozygosity (IHartl and Clarklll997l ). Then the number of sectors will be given by the result for 
the infinite alleles model, Eq. (fT9l) . multiplied by H. 



3.1.3. Finite front size - probability of fixation 

Our previous results were derived under the assumption that the frontier is very large, or rather, 
that sectors are too small to "notice" that the front size is actually limited. However, when a sector 
grows up to a size comparable to the dimension L of the front, we have to account for the possibility 
that the sector could take over the entire front and reach fixation. In this section, we determine 
the probability of fixation assuming periodic boundary conditions, as would be appropriate for a 
radial inoculation, or for microorganisms growing from a linear inoculation around the waist of a 
cylinder. 

As a first step towards this goal, let us determine in the infinite alleles model the probability 
u(r\zo,ro) that a sector has swept to fixation, Z = L, before the frontier has reached position 
r, given that the sector was a size Z(tq) = zo at initial "time" ro- Here, L = 7rd, for a linear 
cylinder of diameter d, and L = 2tt for a radial inoculation. To this end, it is convenient to 
derive from the diffusion equation Eq. (fTTl) an equivalent differential equation for u(r\zo,ro). This 
can be d one, in analogy with derivations of the Kolmogorov backward equation in population 



genetics (IHartl and Clarklll997l ). by writing 



u(r\zo, r ) = J dz u(r\z', r + e)F{z r + e\z , r ) , (20) 

which follows from the Markov property of the problem, see Fig. El Upon expanding the right hand 
side in e and using Eq. (fTTl) . we obtain 



u(r\zo,r ) = J ' dz' ' [u(r\z' \r ) + ed ro u(r\z' \r )] (21) 
x [F(z',r \zo,r ) + 62D(r )^F(z / , r |z , r ) + 0(e 2 )] 
After noting that F(z, ro|^o, ro) = 6(z — zq) and integrating by parts, we find 

d ro u(r\z ,r ) = -2D(r )d^ u(r\z 0: ro) (22) 
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aP(z/a) 




Fig. 6. — Sector size distribution predicted by the annihilating random walker model. The horizon- 
tal axis is the ratio of sector-size z and the positional standard deviation a(r,ro) accumulated by 
a random walk between times tq and r, as defined in Eq. (fT5l) . Full lines represent exact solutions 
with the parameter L/a increasing in steps of 0.5 in order of increasing peak height, as obtained 
from Eq. (|39l) . For L/a > 4, the exact solution is hardly distinguishable from the asymptotic result 
Eq. (fT71) (dashed line), which is independent of L/a. The area under these curves represents the 
probability that a randomly chosen sector has not (yet) reached fixation. For small L/a <C 1, the 
area approaches because the fixation probability of a sector approaches 1. 




Fig. 7. — The graph describes how the probability of fixation depends on the amount of domain 
wall wandering. The function f(a/L) represents the probability that fixation is reached at a front 
of size L given that the accumulated variance of a single domain boundary is a 2 (r, ro), as defined 
in Eq. (US). 
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to order 0{e). Note that all derivatives in Eq. (l22l) act on the coordinates ZQ,ro characterizing the 
initial conditions. 

We seek a solution to Eq. (f22l) subject to the "final" condition u(r\zo,r) = (no fixation if the 
frontier does not advance), and two boundary conditions, 

u(r\0, r ) = u(r\L, r ) = 1 . (23) 

The first boundary condition accounts for the annihilation of a sector as zo — > 0. The second 
guarantees that the fixation probability is 1 for all r > ro if the sector already spans the entire 
front at initial position r = ro- 

Note that the function zq/L satisfies both Eq. ([221) and the boundary conditions Eq. (j23l) , 
which motivates the following ansatz, 

oo 

Au(r\z 0l r ) = u(r\z , r ) - z /L = ^ a n (r )W n (z ) , (24) 

n=l 

with coefficients 

a n (ro) = / dz W n (zo)Au(r\zo,ro) , (25) 
Jo 

given in terms of a complete orthonormal basis set of sine functions 

/2 nn 
- sm(q n z) q n = — . (26) 

The function Au(r\zo, ro) represents the difference between the time-dependent solution, u(r\zo, ro), 
and the linear steady state solution, zo/L. The expansion of Au(r\zo,ro) in terms of sine-modes 
guarantees that the boundary conditions Eq. (f251) are satisfied. 

Inserting the ansatz Eq. ([241) into Eq. (|22l) gives the evolution of the mode amplitudes a n (ro) 
with ro, 

d ro a n (ro) = 2q 2 n D(ro)a n (r ) . (27) 
Upon integrating from r to ro, we obtain 

a n (r ) = a n (r)e"^ 2 /2 ; (28) 

with a = a(r, ro) being the standard deviation defined in Eq. (fT5l) . The pre-factor a n {r) is deter- 
mined by imposing the final condition u{r\zo, r) — on Eq. (l25l) . 



a>n(r) 



[ dzoW n (zo) Au(r\zo,r) (29) 
Jo 

J dz Q sin (^n^o)^ (30) 
L \ — 1 5 otherwise 
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Hence, the solution for the fixation probability of a sector of initial size zq reads 

u(r|*d, ro) = ^ + E sin (q n z ) . (32) 

L n=l ^ 

Clearly, the absolute fixation probability goes to zero as the initial size zo of the sector decreases, 
because small sectors almost always get trapped by the absorbing boundary condition on the left 
side of Fig. 

In the infinite alleles model, we have an initially very large number L/zq of very tiny sectors, 
zo — > 0, that all compete for taking over the colonization front. One and only one sector can reach 
fixation. In order to study the time dependent fixation probability of this successful sector we 
condition on ultimate survival. We ask, what is the probability that fixation is complete by the 
"time" r provided that the conisdered sector reaches fixation? In the infinite alleles model, this 
probability follows from Eq. (|32l) after a normalization, 



/(rr ) = hm— , r (33) 

z ^o u{a —> oo|zo, Tq) 



= l + 2^(-l) n e-^ fj2 



/2 



n=l 



= tf 4 [0,exp(aV/(2L 2 ))] , (34) 

where $ n [z,q] are the Jacobi theta functions. The function f(r\ro) represents the probability that 
a colonization experiment running from time ro until time r reaches fixation in the infinite alleles 
model. The denominator in Eq. ([331) is a normalizing factor that ensures f(r\ro) —> 1 as a(r) — > oo. 
The dependence of f(r\ro) on a/L is shown in Fig. [71 

At linear fronts, the accumulated standard deviation a = ax oc y/Ar increases without bound, 
so that fixation is inevitable at long times. For this case, (again restricting our attention to the 
infinite alleles model) we can determine the average "time" or frontier distance to fixation, (Arj^), 
by an integral over /(r|ro), 

POD 

(Ar fix ) = / dr (r - r )d r f(r\r ) (35) 

Jr 

oo 

dr [1 - /(r|r )] (36) 
= 5 (linear inoculations) (37) 

where we used Eq. (OQ) for the variance o\ and Eq. f[34l) to evaluate the integral in the second line, 
and the boundary condition lim r ^ 00 f(r\ro) — > 1. Note that the pre-factor in Eq. (j37l) holds for 
periodic boundary condi tions only, bu t can be evaluated for other boundary conditions along the 



ro 

2 



same lines, see e.g. Ref . ([Rednerl 120071 ) . Our derivation assumes the infinite alleles scenario, and 



we have not yet found an exact result for a finite number of colors. We expect, however, that the 
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fixation time is linearly dependent on the heterozygosity H of the initial population, at least for 
small heterozygosities. Our approximate argument is as follows. Suppose the initial population 
is an unbalanced binary mixture with H <C 1. Then, with high probability (1 — H) the majority 
population will take over in a short time of order (HL) 2 /D. On the other hand, with a low 
probability H the colony will be taken over by the minority population, and this will take a time 
comparable to the fixation time L 2 /D in the infinite alleles model. As a consequence, these rare 
events of the fixation of the minority population dominate the avarage fixation time, which should 
therefore be of order HL 2 /D. 

Finally, we also generalize the sector size distribution Eq. (fTTI) to finite-sized frontiers. To this 
end, we solve Eq. (fTTI) under an additional absorbing boundary condition, 

F(z,r\z ,r )\ z=L = (38) 

which amounts to disregarding a sector once it reaches fixation. In Supplementary Text 1, we show 
that the solution to this problem can again be found in terms of the sine modes introduced in 
Eq. (j26l) . specifically 

F(z, r\z , r ) = — ^ sm(q n x) sin(g n z )e"^ a2/2 . (39) 

71=1 

As in the unbounded case, the sector size distribution is found by letting zq —> 0, 

00 

P(z, r\r ) oc ^ q n L sin(g n x)e-^ 2 / 2 . (40) 

71=1 

In Fig. [6J we plot this sector size distribution normalized by the probability that a sector has not 
(yet) reached fixation. The results are indistinguishable from the asymptotic result Eq. (fTTI) for 
L/a < 4. 



3.2. Natural selection 

The above model can be extended to describe the effect of natural selection on mutations 
occurring during a range expansion, such as the one documented in Fig. 03 For simplicity, we 
remove complications associated with inflation by focusing on linear inoculations. In the previous 
section on neutral evolution, we ensured that each sector had the same chance of survival by 
requiring that the expected change in sector size vanishes, (AX) = 0, see Eq. (pQ). This neutrality 
assumption is no longer valid if a sector is generated by mutants that have a different fitness than 
the wild type. To describe the fate of those non-neutral sectors, it is natural to assume that a 
selective advantage biases the sector growth. A sector harboring beneficial mutations will tend to 
increase its size, whereas a deleterious mutation will decrease the size. Upon setting Ar = r — ro 
as before, we thus generalize Eq. (pQ) to 

(AX) = 2m ± Ar , (AX 2 ) = 4D x Ar . (41) 
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Here, we have assumed that selection does not alter the diffusion constant Dx, which seems reason- 
able, at least in the case of weak selection. The new parameter 2m± describes the increase in the 
mean sector size due to selection and has units of a slope, length/length. The factor of 2 indicates 
that each of the two boundaries exhibits an average lateral drift of raj_, adding up to a total sector 
growth rate of 2m±. In terms of measurable units, the bias parameter m± is given by the ratio of 
the velocity v_\_ of a sector growth at right angles to the direction in which the front is advancing 
and the velocity v of the wild type range expansion, 



m_L = — = tan($/2) . 



(42) 



The second equality describes the relation between m± and the opening angle $ of the sector at long 
times, which can be perceived from the sketch in Fig. [91 For weak selection, $ will be small, and we 
can think of m± as being just half the (asymptotic) opening angle of the sector. This opening angle 
depends on the relation between the fitness effect of the mutation and the demographic expansion 
process. 

Although the relation between $ (or m±) and selective advantage is complicated in general, 
it can be determined in two simple cases, both of which should be realizable for populations of 
microorganisms on a Petri dish, and for other range expansions as well. We assume that the only 
phenotypic effect of a beneficial mutation is an increase in the expansion velocity v — ► v* — f(s)u : 
where f(s) > 1. Here, s is the usual selective advantage, defined as the ratio of growth rates a* and 
a of a population of organisms at the frontier, a* /a = 1 + 5. In principle, a (or a*) could be measured 
directly by monitoring cell divis ions under a microscope, as was done in the sectoring experiments 
of Ref. (jHallatschek et al.1 120071 ). However, the relation between the more easily measured front 
growth velocities and s is known when population number fluctuations (i.e., effects of genetic drift) 
at the front are weak . One then expects the range expansions to be described by Fisher population 
waves (jMurrayl 12004 ) , for which v = 2 V Da and v* 



2\fD*a* , where D and D* are population 
diffusion constants at the frontier. If the beneficial mutation does not change the diffusion constant, 
we then have 

/(*) = V^T^ . (43) 

Wakita et al. have in fact measured a squa re root dependence of the growth velocity on the nutrient 
concentration in plates of bacillus subtilis (IWakita et al.lll994l ). consistent with Eq. (f43l) above when 
the growth rate is proportional to the nutrient conc entration. Results a re also known in the limit 
when genetic drift dominates growth at the frontier (jDoering et al.ll2003l ). In this limit, the growth 
velocities are proportional to doubling rates at the frontier. Assuming all other quantities are the 
same for the mutant and wild type, we then have 

f(s) = 1 + 8. (44) 

See Supplementary Text 2 for a more detailed discussion. 

More generally, we expect that f(s = 0) = 1, and a Taylor series expansion about of the form 

f(s) = l + cs + ... (45) 
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where c > and s are small. Under the conditions described above, we have c = 1/2 and c = 1 for 
weak and strong genetic drift respectively. 

Now suppose that a mutation arises (or a mutant cell is inserted by hand) close to the linear 
front of a wild type population in a Petri dish, and is able to overcome the critical initial phase, 
where selection is weak compared to genetic drift. For long times, the mutant sub-population will 
then form a sector that grows by the factor f(s) faster into the unoccupied space than the wild 
type, as discussed above. Assume that, after initial transients, a time-independent sector angle 
forms. Again appealing to Fig. [91 we see that the mutant sector asymptotically approaches an 
opening angle $ that satisfies 



Thus, for this "geometric" model of the opening angle the drift parameter m_\_ = v_\_/v in Eq. f|42l) 
takes the form 



where the approximations above assume s < 1. 

It seems likely that this simple phenomenological model applies to the opening angles created 
by mutant strains of microorganisms on a Petri dish. Note that the square root dependence in Eq. 
(51) suggests that the sector angle could be a quite sensitive measure of weak selective differences. 
Other functional relations between and s are possible, however. As discussed in Appendix B, when 
number fluctuations are strong compared to the selective advantage at the front (strong noise limit), 
a linear relation arises in two-dimensional stepping stone models. These models allow some inter- 
diffusion of mutant and wild type strains after the population wave has passed by. A linear relation 
can also arise for the stochastic Fisher genetic waves generated associated with weakly deleterious 
mutations. 

It remains to be seen experimentally which model most accurately describe the function m±(s) 
in a given situation. However, our parameterization of beneficial mutations in terms of the phe- 
nomenological parameter m± is quite generally applicable to any growth scheme leading to a bias 
m± in the random walk of domain boundaries. In the absence of a detailed microscopic understand- 
ing, one can always choose to parameterize beneficial and deleterious mutations directly in terms 
of m± itself. Although the arguments above focus on m± > 0, it is easy to see that deleterious 
mutations are described by m± < 0. In the following, we analyze the evolutionary consequences of 
this bias for both positive and negative values of m±. 



cos($/2) = — = 



1 — CS . 



(46) 




(47) 



3.2.1. Beneficial mutations, m± > 



Isolated beneficial mutations are often lost due to genetic drift. However, because a sector of 
beneficial mutation tends to increase in size, m± > 0, its survival probability is increased compared 
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Fig. 8. — Graphical representation of Eq. (f20l) for the probability u(r\zo,ro) a sector fixes at any 
frontier position up to and including r, given that the width Z of the sector was zo at earlier frontier 
position ro- Fixation occurs when the fluctuating width (black line) of a sector matches the system 
size L. 
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Fig. 9. — A simple model of a beneficial sector within a planar wave front in the long-time limit. We 
assume that a beneficial mutation arises at front position ro and is able to overcome the short-time 
genetic drift at the front causing wobble in the domain boundaries. At long times the beneficial 
sub-population segregates and forms its own sector growing at radial velocity u*, which is assumed 
to be larger than wild- type growth speed v. Requiring that the kinks at the interface, where the 
wild type and mutant population meet, are reached by both populations at equal times, leads to 
cos(<3>/2) = v/v*. 
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to the neutral case. Here, we determine the fixation probability of a sector of beneficial mutations, 
and derive how frequently those sectors appear, given a certain rate of beneficial mutations. 

Similar to the unbiased case, we will determine the fixation probability from a diffusion equation 
that describes the statistical properties of the size X of a sector. Let u(r\xo,ro) be the probability 
that a sector of beneficial mutations reaches fixation up to and including front position r given that 
it had size X(ro) = xo at initial position tq. Using analogous arguments to those in Sec. 13.1.31 it 
can be shown that the assumptions in Eq. (I4T1) lead to a biased diffusion equation for the fixation 
probability u{r\ XQ,ro), which reads 

d ro u(r\x , r ) = -2D x d% Q u(r\xo, r ) - 2m±d Xo u(r\xo, r ) • (48) 

Apart from the new drift term oc m± this equation is identical to the unbiased case Eq. (l22l) with 
the constant diffusion parameter appropriate to linear inoculations. We again impose the boundary 
conditions of Eqs. ([231) accounting for annihilation at x = and fixation at x = L, with the periodic 
boundary conditions appropriate to growth along a cylinder. 

Because the diffusion constant and drift parameter are both independent of the time-like 
frontier position variable in the linear inoculation, we seek a ro-independent solution of Eq. (|22l) . 
After setting the right-hand side to zero and integrating twice, we find the steady state solution in 
the limit of long times, r(t) — > oc, 

I _ e -m ± x /D x 

U||(oo|£ ,ro) = 1 _ e . m±L/Dx > (49) 

which is the ultimate survival probability of a beneficial mutation at a linear front 0. 

The exponential dependence of Eq. (JJ9|) on xq implies that a sector almost certainly overcomes 
stochastic loss when it reaches a size larger than an "establishment length" / = Dx/m±. If a sector 
is much smaller than this characteristic length, xq <C /, the survival probability takes the simple 
form un « m±xo/Dx provided that L/l 1. This result can be used to relate the frequency at 
which beneficial mutations become established in the form of sectors to the beneficial mutation rate 
jlbj which has units of an inverse time. To this end, we assume that all beneficial mutations confer 
the same selective advantage raj_, and consider the evolutionary dynamics during short effective 
"time" increments Ar, in which genetic drift is stronger than selection. According to the neutral 
results of Sec. 13.11 only a number N(Ar) of front lineages evade stochastic loss during a time 
increment Ar. Equivalently we can say that there is at any time a set of 7V(Ar) individuals whose 
descendants will be present after the next time increment Ar. Among this population of founders, 
beneficial mutations occur at rate /i&iV(Ar) 5 where fi^ = jib/v is the beneficial mutation rate in 
units of an inverse length. These mutations will ultimately survive at long times with probability 



4 Eq. (|49J is formally sim ilar to Kimura's fixation probability of a beneficial mutation in a well-mixed popula- 
tion ([Crow and Kimurall 19701 ) when we identify the product of population size and selection coefficient with m±L/2Dx 
and the frequency of the beneficial allele with x/L. 



-23- 



m±xo/Dx, where xq = L/N(Ar) is the average size of a sector after time Ar. Multiplying these 
factors together yields the rate at which beneficial mutations become established 

Beneficial mutations establishment rate = ji^Lm^j Dx • (50) 

Note that the Ar-dependent sector number N(Ar) drops out of the final result H. Note also that the 
establishment rate is an inverse length in our choice of units, which has the following interpretation: 
It measures the number of beneficial mutations appearing per unit length of progression of the front. 
As discussed below, in situations where beneficial mutations are of the type proposed in Fig. [9J 
Eqs. (|46l) and ([501) can form the basis of an experimental technique to measure the relative fitness 
advantage of beneficial mutations (as embodied in the parameter m±) as well as beneficial mutation 
rates in microbial colonies. 



3.2.2. Deleterious mutations, m± < 



Given the importance of chance effects during population expansions, one may wonder to 
what e xtent deleterious m utations can prevail at expanding frontiers. Indeed, a recent simulation 
study (ITravis et al.1 120071 ) has observed deleterious mutations that are swept to high frequencies 
by population waves. It is easy to see how these "gene surfing" events of deleterious mutations 
may come about in one spatial dimension. The first step is a matter of chance. A deleterious 
mutation needs to arise close to the wave front and become frequent there, despite natural selection. 
Such surfing attempts are promoted by the stron g genetic drift at expanding fr ontiers, which are 
characterized by small effective population sizes (IHallatschek and Nelson! 120081 ) . Once the front 
has been taken over by a deleterious mutation, the population wave advances into empty territory 
at a somewhat reduced velocity v* < v due to the reduced fitness of the mutants heading the 
expansion. The wild type population, on the other hand, is stuck in the bulk of the population, 
but nevertheless advances by displacing the less fit muta nts. The ensuing genetic wave of advance 
may be described by a Fisher genetic wave (lFisherlll937l ) with velocity v g (see Fig. [TO]). Although 
the wild type population tries to catch up with the front via this genetic wave, this will actually 
never happen if the genetic wave is slower than the population wave of the mutants, v g < v*. Thus, 
deleterious mutations may take over permanently if v* > v g . See Supplementary Text 2 for an 
example of a one dimensional model that exhibits this behavior over a broad parameter range. 

In two spatial dimensions, isolated deleterious mutations cannot permanently prevail at ex- 
panding frontiers. Although a sector harboring deleterious mutations may arise by chance effects 
as well, this sector is ultimately doomed to extinction because it has a lower expansion velocity 
than the surrounding wild type sectors. As a consequence, it is just a matter of time until the 



5 Note that our argument assumes that the width of sector boundaries is smaller than the establishment length /. 
Otherwise, it is not possible to describe the sectoring dynamics as a random walk of boundaries when they are closer 
than /. 
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deleterious mutants are literally overtaken by the wild type population, see Fig. [2] and the time 
sequence illustrated in Fig. [TTk . Nevertheless, deleterious mutations that temporarily surf on a 
population wave - until they finally "fall off" the wave front - pose a potentially serious threat to 
the pioneer population because they achieve much higher frequencies than expected under well- 
mixed conditions. As we demonstrate below, when these deleterious alleles become numerous, they 
can even trigger mutational meltdown in a cooperative manner. 

We seek to quantify the frequency of surfing deleterious mutations at expanding frontiers within 
the diffusion approximations used in this paper. As before, let F(x, r|xo, ro) be the probability that 
a deleterious sector has size X(r) = x at frontier position r, given that it had size X(tq) — xq at 
earlier "time" tq. The biased diffusion equation for this distribution function F(x,r|xo,ro) reads 



This equation generalizes the unbiased Eq. (fTTI) . formulated backwards in time, and has the same 
form as Eq. (1481) for the fixation probability u(r\xo,ro) of a beneficial sector. The only difference 
is that the drift term has the opposite sign describing the shrinking of deleterious sectors. 

We will now use equation Eq. (|5T1) to study the statistical properties of the area A{x$) spanned 
by a deleterious sector of initial size xo (see Fig. [TTb). When combined with the mutation rate, 
this quantity turns out to control the steady state frequency of deleterious mutations at an advanc- 
ing population front. On the "space-time" plot of two random walkers, where r is the time-like 
coordinate in the growth direction, A{x$) simply represents the area enclosed by the two colliding 
world lines, given they are initially separated by a distance x$. The average area (A(xq)) can be 
expressed as an integral of F(x,r|xo,ro) over its final coordinates, 



To obtain a differential equation for (A(xq)), we multiply Eq. (|5T1) by x and integrate over x and r, 



d ro F(x, r\x , r ) = -2D x dl Q F(x, r|x , r ) + 2 \m±\ d XQ F(x, r\x , r ) . 



(51) 




(52) 



roc /•( 

/ dr / 

Jr JO 



■oo 



dxxd ro F(x,r\xo 1 r ) = -2D x d 2 XQ (A(x Q )) + 2 \m±\ d XQ (A(x )) . 



(53) 



WT 




Fig. 10. — A deleterious mutation can permanently take over a population front in one spatial 
dimension if the speed of the genetic wave v g of the wild type (WT) invading the mutants (MT) is 
lower than the speed v* of the advancing population front of mutants. 
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Because the left hand side represents the total derivative d ro (A) up to a boundary term evaluated 
at r = ro, we have 



d ro (A(x )) + / dxxF(x,r \x ,r ) = -2D x d 2 (A(x Q )) + 2\m±\d XQ (A(x )) 



(54) 



Upon noting that F(x, ro|xo> ro) = 6(x — xq) and that the expected mutant area (A) = (A(xq)) is 
ro-independent in linear inoculations, we obtain a simple ordinary differential equation for (A(xq)), 



2D x d 2 XQ (A(x )) - 2 \m ± \ d Xo (A(x )) + x = . 



This equation is solved by 



(A(x )) = 



Xn 



4 |m_L 



DxXq 

2m\ 



(55) 



(56) 



The quadratic contribution to the area, which dominates for xq ^> I = Dx/\m±\, is just the 
deterministic expectation for the area of the mutants neglecting genetic drift: With no diffusion 
of its boundaries, the deleterious sector should shrink laterally at a rate 2|m^| and thus collapse 
when the front has advanced by r c — xo/(2\m±\). The deterministic sector area should thus equal 
the area XQ/4|m^| of an isosceles triangle with height r c and base xq. Note that the characteristic 
length I = Dx/\m±\ now plays the role of a "disestablishment length". 

The linear part in Eq. (l56l) is due to stochastic fluctuations and dominates on small scales, 
for £o <C I. As we now show, this stochastic part determines the mutational load of an expanding 
population. As in Sec. 13.2. 1\ we note that stochastic fluctuations dominate over selection on a 
small "time" scale Ar. There are at any time N(Ar) founders - mutants or wild type - that give 
rise to sectors of average size xq = L/N(Ar) at a small time Ar in the future. Among these 
N(Ar) individuals deleterious mutations occur at a rate A/"(Ar)/^, with ji^ — jld/ v representing 
the effective deleterious mutation rate per unit length of frontier growth. This process produces 
deleterious mutants at a rate given by 



fi d N(Ar) (A[L/N(Ar)}) = 7 L 



(57) 



t ■ 





6 












Fig. 11. — a) A deleterious mutation (red) arising close to the population front (black lines) achieves 
high frequency by temporarily "surfing" at the front. Figure b) defines the quantity A{x$), which 
measures the expected area generated by a mutant sector if it has initial size xq. 
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For Ar — > 0, the short time neutrality assumption becomes exact, the factors of N(Ar) cancel and 
we obtain 

DxHd /rox 

7 = i^T (58) 

The restrictions of the footnote after Eq. ([501) apply here as well, and we have assumed for simplicity 
that all deleterious mutations confer the same selective disadvantage — |raj_|. 

The parameter 7 characterizes the mutational load. If 7 <C 1, it represents the fraction of 
the individuals in the newly colonized regions that carry the deleterious mutation. When the 
parameter 7 becomes of order one, mutant sectors become so numerous that they start to collide. 
These collisions are not captured by our simple theory, so 7 can no longer be interpreted literally as 
the fraction of mutants for 7 = O(l). Still, 7 defined by Eq. ([581) can be viewed as a dimensionless 
control parameter that determines whether the pioneer population remains close to wild type or 
the mutational load becomes so strong that the average fitness deteriorates dramatically. 

In fact, simulations for the simple model described in Supplementary Text 3 do indeed reveal 
a critical value 7 C = 0.32 ± 0.02 that separates two regimes, see Fig. [12l For 7 < j C: the average 
fraction of wild type surviving in a pioneer population stabilizes at a finite steady state value, 
whereas for 7 > 7 C the pioneer population keeps accumulating deleterious mutations until the wild 
type is lost entirely. This apparent continuous phase trans ition is the one dimensional spatial analog 
of the error threshold in well-mixed populations (jM.lll97ll ). Note that, depending on the functional 
relation between m± and the selective disadvantage s < 0, the threshold value 7 C for expanding 
populations can be very different from the well-known prediction <^™- m - — = 1 for well-mixed 

populations. For the "geometric model" outlined at the beginning of Sec. 13.21 we have m± oc yfs 
and thus 7 C oc s -1 similar to a well-mixed population. However, for a two-dimensional stepping 
stone model in the strong noise limit, the relationship between m± and selective disadvantage s 
is linear (cf. Supplementary Text 2). In this case, the genetic load 7 oc s~ 2 could be strongly 
increased due to weakly deleterious alleles. 

A simple argument suggests that the collision of mutant domains will indeed increase the 
expected fraction of mutants. Consider the collision of two mutant domains of sizes x\ and £2, 
characterized by the same selective disadvantage \m_\_\. Then the newly combined domain has the 
increased size x\ + X2. The expected area of the combined region as given by Eq. ([561) is larger 
than the two separated domains by an amount 

(A ex (x u x 2 )) = (A( Xl + x 2 )) - (A( Xl )) - (A(x 2 )) = . (59) 

This excess mutant area generated by the collision of domains furnishes a cooperative interaction 
associated with the mutational meltdown, as revealed in Fig. [12] as 7 — > j c . 
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Fig. 12. — Simulation results for the fraction pwr of the wild type population for various values 
of the drift parameter |raj_|, characterizing the selective advantage of the wild type population, 
as a function of the dimensionless control parameter 7 = /i^Z?x/(2m^). The simulations were 
carried out for small values of m± because our theoretical treatment applies to the limit of weak 
selection. To limit computational intensity, we restricted the simulations to m± > 0.01. Upper 
figure: With decreasing raj_, Pwr{l) seems to approach a universal function. The simulations 
confirm our analytical result that 1 — pwr ~ 7 as 7 — » 0, see inset. On the other hand, they also 
show that pwr — ► above a critical threshold. The lowest probed value of m± = 0.01 indicates 
7 C = 0.32 ± 0.02. The simulation algorithm is described in Supplementary Text 3. 
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4. Conclusions and discussion 



We have studied the impact of a continuous range expansion on the evolutionary dynamics of 
spatially structured populations. Our analysis focuses on the enhanced genetic drift at advancing 
frontiers, which leads to strong genetic differe ntiation. Most directly, our results apply to recent 
experiments on expanding microbial colonies (IHallatschek et al.l 120071 ). However, we believe that 
the present analysis is of more general significance for populations that grow continuously and 
isotropically in two dimensions. The advancing frontier of those populations also has a thin layer of 
active pioneers, whose width depends on demographic parameters of the range expansion, such as 
the growth, migration and turnover rates. As the population advances, the pioneer population in 
this quasi-on e dimensional front region is continuously re-sampled, as in a one-dimensional stepping 
stone model (IKimura and Weisslll964l ). for which it is known that allelic segregation occurs locally 
beyond a crossover time[§. Therefore, we believe that the coarsening process, and our model for it 
in terms of sharp domain boundaries, might be generally important for homogeneously expanding 
populations. 

Within this model for population dynamics and genetic segregation of a continuous two- 
dimensional range expansion out of some prescribed initial habitat (linear or circular), we first 
studied the neutral genealogies of single loci without mutations. We found that only a fraction 
of founder cells are able to propagate their genes with the advancing front. Due to population 
turnover in the thin band of pioneers, the number of lineages gradually decreases as the expan- 
sion progresses. If a linear front advances by a distance Ar, the number of survivors decays like 
1/y/Ar. One sector will dominate after an average fixation time oc L 2 , where L is the linear di- 
mension of the front. In the case of circular colonies, we find instead that a finite number of sectors 
survives, due to the geometric expansion of the perimeter which opposes genetic drift. For our 
moel, the expected number of these surviving sectors is proportional to the inverse square root of 
the initial radius of the colony. Both results assume that domain boundaries carry out a diffu- 
sive random walk. Growt h models with surface roughness are expected to show an a nomalously 
vigorous wall wandering ( Saito and Muller-Krumbhaar 19951 : Hallatschek et al. 2007 ) where the 
variance in separation X(r) between two walls grows as (Ar) 2 ^ with an exponent 2( > 1 as the 
front advances by a distance Ar. Our analysis of annihilating random walkers could be extended 
to these cases by assuming a diffusion "constant" that decreasing with the separation between the 



6 In a continuous one- dimensional stepping stone model, well-defined boundaries appear after the characteristic 
time T — Ci d DTg, where D is the usual diffusivity of the individuals (with dimensions cm 2 /sec), r g is the generation 
time characterizing the strength of population turnover, and c±d is the one-dimensional population density with units 
of an inverse length ( Barton et al.ll2002 : Hallatschek and Korolevll2008 ). Considering the advancing edge of the colony 
as an effectively one-dimensional habitat of width v/a depending on the expansion velocity v and effective growth rate 



a leads to the estimate c\a ~ c^dvja. Thus, we expect segregation to occur on times larger than T « Dc 2 



Domains will then occur on length scales larger than w 
width of domain boundaries. 



>DT « Dc2dVT g /a, where w represents a characteristic 



-29- 



two random walkers 0. On the scaling level, such an analysis would yield an average sector num- 
ber for linear inoculations increasing like (Ar) - ^ instead of (Ar) -1 / 2 . For radial inoculation, the 
number of surviving sectors 7V(oo,ro) (see Eq. (j7j)) now scales with the initial radius according to 
7V(oo, ro ) ~ r o~^- We speculate that the effect of inflation balancing genetic drift could be relevant 
for explaining the surprisingly large levels of genetic diversity among some invasive species, that 
have been introduced locally into a new and favorable habitats (e.g. rabbits in Australia). Ac- 
cording to our model for radial expansions, those species are expected to preserve a considerable 
amount of their initial diversity during the habitat expansion. 

Next, we considered the evolutionary dynamics of beneficial mutations arising at frontiers. To 
this end, we modified our model by adding a deterministic bias to the sectoring dynamics that 
tends to increase the size of sectors. This modification applies most directly to beneficial mutations 
whose only effect is to increase the growth rate of individuals v — > v* = f(s)v. We found that, at a 
linear frontier, a successful beneficial mutation "emits" a sector that has an opening angle $ given 
by Eq. (H6l) . If the change in velocity due to the mutations is weak, s <C 1, then the sector angle 
is approximately given by $ « 2a/2cs, where c — 1/2 or c — 1 depending on wether genetic drift is 
weak or strong at the front, respectively. The square root dependence on the selection coefficient 
indicates that sector angles are a sensitive measure of fitness differences. Thus, measuring sector- 
angles could be a useful tool to decipher the distribution of fitness effects of beneficial mutations, 
which can be difficult (or tedious) to measure by liquid culture techniques. The use of advancing 
population waves in evolutionary studies has been demonstrated ea rlier in a one-dimensional st udy 



of the effect of beneficial mutations on replicating RNA molecules (IMcCaskill and Bauerlll993l ). In 
this study, beneficial mutations were identified by a spontaneous increase in expansion velocity of a 
Fisher population wave of proliferating RNA molecules travelling along a one-dimensional tube. In 
the two-dimensional extension of these studies, we expect a marked increase in resolution because 
we expect sector angles to be much more sensitive to selective changes than spreading velocities. 

Although our parameterization of beneficial mutations in terms of the phenomenological pa- 
rameter m_\_ is quite generally applicable, the relation between bias m± and fitness difference 
remains to be assessed experimentally. The example of a beneficial mutation in Fig. [3] indeed seems 
to approach an asymptotic sector angle. The sector shape deviates initially somewhat from the 
simple triangular sector geometry used in our purely geometric model, however. This transient fea- 
ture suggests that a number of other factors could change the relation bet ween bias m± for small 



sectors and the fitness difference in microbial populations (IRon et al .1120031 ). such as a line tension 
penalizing front deformations. A related possibility is an effective attraction between oppositely 
directed Fisher genetic waves that slows down their separation at early times. 

Finally, we estimated the mutational load in an expanding population. Deleterious mutations 
were found to proliferate at expanding frontiers rather differently than in well-mixed populations. 



7 At a linear front, for instance, the diffusion constant Dx(Ar) ~ (Ar) 2 describes a super-diffusive random 
walk with exponent f . 
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Given that deleterious mutations with a selective disadvantage \m±\ occur at a small rate /i^, 
we determined the fraction 7 of mutants at a linear wave front as 7 = Dxfid/ \2vrrij_). Here, v 
is the expansion velocity of the population wave, Dx is the wall diffusion constant (in units of 
length 2 /length) and m± represents the "speed" or slope by which a single deleterious mutation 
is squeezed out of the population frontE]. Furthermore, we presented numerical evidence that 
populations suffer from genetic meltdown as 7 — > 7 C « 0.32 ± 0.02. This error thres hold should 



be compared with the expectation 7^. m . = jjid/ s f° r well-mixed populations (lM.lll97ll ). where s is 
the relative fitness detriment of the deleterious mutations. Depending on the product j/j w .m. = 
Dxs/(2vm 2 L ), both predictions can be very different. If wall wandering represented by Dx is 
strong, or the bias m± for a given fitness difference is weak (for instance m± oc s as in the weak 
selection limit of the stepping stone model, cf. Supplementary Text 2, then the genetic load during 
a range expansion becomes substantially larger than in the well-mixed case. Thus, deleterious 
mutations could accumulate quite strongly during times of habitat expansions, thereby setting 
tight constraints on the dynamics of range expansions. Genetic load could be serious threat during 
general species invasions. This effect might extend to past range expansions of humans, as the 
mutatio nal load inside Africa wa s found to be significantly lower than in the more recently colonized 
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A. Diffusion equation for periodic boundary conditions 

In the main text, we gave a simple argument how to solve the diffusion equation (fTTl) for 
annihilating random walks with absorbing boundary conditions at z = 0. However, this solution 
is only strictly valid in unbounded space. Here, we derive the slightly more complicated solution 
Eq. (1391) valid for finite systems. This solution satisfies another absorbing boundary condition, 



8 |m^| is approximately given by half the "closing" angle of a sector harboring deleterious mutations. 
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Eq. (l38lh which accounts for the fixation when the sector reaches the size of the system L. 
We make the following ansatz 

oo 

= 5> n (r)W n (z) (Al) 

n=l 

a n {r) = [ dzW n (z)F(z,r\z ,r ) , (A2) 
Jo 

where the sine mode W n (z) and wave numbers q n were defined in Eq. ([26]) . This expansion in terms 
of sine-modes on the right-hand-side guarantees both boundary conditions. 

Observe that the ansatz Eq. fjAll) solves Eq. ([TTD provided that the mode amplitudes a n (r) 

obey 

d t a n {r) = -2q 2 n D(r)a n (r) . (A3) 

Integrating this from ro to r gives 

a n (r) = a n (r )e-^ 2 / 2 (A4) 

where a is the standard deviation in the z coordinate, as defined in Eq. (fT5l) . The pre-factor is 
determined from the initial condition F(z, ro\zo, ro) = S(z — zo) : 

a n (r) = [ dzW n (z)F(z 1 r\z 0l r ) (A5) 
Jo 

|- sm(q n z ) . (A6) 
Hence, the solution for the probability distribution reads 

F(z, r\z , r ) = j ^ sm(q n z) sm(q n z )e- q ^ a2 /2 , (A7) 



which is Eq. (1551). 



Next, we would like normalize the distribution F(z, r\zo, ro) by the total probability ^(r|zo 5 ^0) 
that a sector of initial size zq has not collapsed. This probability can be written as a sum of two 
terms, 

u(r\z 0l ro) = ^(r|z ,^o)+/ dz F(z, r\z , r ) (A8) 

Jo 

The first term on the right hand side in Eq. ()A8j) is the probability that the sector has reached 
fixation, which was evaluated in Eq. (f32l) , and the second term is the total probability that a sector 
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has neither collapsed nor reached fixation. In the limit zo — ► 0, we can now construct the normalized 
distribution 

P(z r\m) - lim F(z ' r|z °' ro) - Si 2q n sm(q n z)e-^ 2 / 2 

P(z,r\r )- hm - ^ [Q> ^ (<Ara/(2L2))] • (A10) 

This limiting distribution is depicted in Fig. [6] for various values of L/a. Note that the results 
are indistinguishable from the asymptotic result Eq. (fT71) for L/a < 4. The area under these 
plots represent by construction the probability that a randomly chosen sector has not yet reached 
fixation. 



B. Genetic versus population waves 

In this Appendix, we first exhibit an explicit model for the surfing of a deleter ious gene in one 



dimension. Unlike the two-dimensional case considered in Ref. iTravis et al.l (120071 ). we focus here 
on one dimension, leading to the situation sketched Fig. [101 Provided we focus on populations that 
stop growing and diffusing once the population has passed by when c + c* = 1 (see below), we 
can view this model as an approximation to the dynamics within the thin layer of actively growing 
pioneers in a two-dimensional range expansion, as in FigOJi) and[T]i). This approximation neglects 
number fluctuations at the front. The time variable t in this Appendix then corresponds to the 
frontier position r used elsewhere in this paper. 

Let c(x,t) be the dimensionless concentration of wild type individuals (growth rate a), and 
c*(x,£) the dimensionless concentration of the mutant strain (growth rate a*). We assume for 
now < a* < a, so that the mutation is deleterious relative to the wild type. We work in the 
strong selection limit, and neglect fluctuations in the number of discrete individuals within the two 
populations, although these can be important under some circumstances. For simplicity, we assume 
identical diffusion constants and a common steady state value or "carrying capacity" of unity in 
rescaled units for these two populations, both separately and when they are mixed together. The 
two populations differ in their growth rates, and in addition the wild type secrets a chemical that 
impedes the growth of the mutant under crowded conditions. A simple set of coupled reaction- 
diffusion equations for the two strains then reads, 

d t c = Dd 2 x c + ac(l-c-c*)-ac*c, (Bl) 
fyc* = Dd^c* + a*c*(l - c* - c) + acc* . 

"Crowded conditions" corresponds to c + c* = 1, and the term —ac*c (with a > 0) in the first 
equation represents the secretion of an inhibitory chemical. Consider first the zero-dimensional 
case of well-mixed, spatially uniform populations, so that c and c* are a function of time only. It 
is straightforward to show that evolution of c(t) and c*(£) is then controlled by three fixed points, 
namely 



(0,0), eigenvalues a Sz a* (unstable) 
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(0, 1), eigenvalues — a & a (hyperbolic) 
(1,0), eigenvalues — a* & — a (stable) 



(B2) 



and that the invariant subspaces are c = 0, c* = 0, c + c* = 1. The equal and opposite coupling 
strengths of the last terms in Eqs. (IB1I) were chosen to insure that c+c* = 1 is an invariant subspace, 
and that we deal with the relatively simple situation that a composite population has the same 
carrying capacity as either population in isolation. The dynamics associated with this well-mixed 
effectively zero-dimensional population is shown in Fig. [131 If a* ~ Q> > 0, both populations will 
grow up rapidly after a small inoculation near (0,0). However, provided < a C 1 , the ultimate 
fate of the population is then a slow drift down the line c + c* « 1 until the wild type dominates 
at the stable fixed point (1,0). 




Fig. 13. — Dynamics associated with Eqs. (IB2I) in zero dimensions, e.g., in a well-mixed overnight 
culture of microbes where spatial diffusion plays no role. 

Now consider what hap pens in one spa tial dimension. If the mutant strain is absent, (c*(x, t) = 
0), standard considerations iMurrayl (|2004l ) show that the remaining equation, 



d t c = Dd 2 x c + ac{l-c) 



(B3) 



supports stable left- and right-moving waves that interpolate between c(x,t) = and c(x,t) = 1 of 
the form 

c(x,t) = f(x±vt) , (B4) 



with the velocity 



v 2VDa 



(B5) 



and width w = yjD/a. On the other hand, if c{x,i) — everywhere, the Fisher wave describing 
the mutant population c*(x,t) has velocity 



(B6) 
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Not surprisingly, the mutant Fisher population wave spreads with a lower velocity than the wild 
type, v* < v. However, consider now the situation when both strains are present and together always 
saturate the carrying capacity of the environment, i.e., we operate along the line c(x, £)+c*(x, t) — 1. 
Equations Eq. (jBll) then collapse to a single equation. If we focus on the dynamics of the wild type 
density c(x,t), we now have 

d t c = Dd 2 x c + ac{l-c) . (B7) 



Equation Eq. (IB7I) describes a Fisher genetic wave of the wild type displacing a saturated population 
of mutants with a velocity 

v g = , (B8) 

determined by a, which represents the selective advantage under crowded conditions. Note that 
the Fisher genetic wave velocity is proportional to the square root of this selective advantage a. 

Figure [10] depicts in one dimension a superposition of a Fisher population wave on the right, 
with a mutant population moving into empty space with velocity u*, and a Fisher genetic wave on 
left, with the wild type displacing the mutant with velocity v g . Successful "surfing" of a deleterious 
mutant ahead of the wild type requires only that the genetic wave not overtake the population 
wave, i.e., v g < v*, which leads to the condition 

a < a* < a (B9) 

for this simplified model. 

As discussed in the main text, the situation is considerably more complicated (and interesting!) 
in two dimensions, and when the genetic drift embodied in particle number fluctuations are taken 
into account. Let us focus now on a beneficial mutation forming a sector like that in Fig. [91 We 
consider the simple situation such that a* w a (comparable Fisher population wave velocities), 
but assume that the wildtype is favored under crowded conditions, a < with |a| ^> \a — a k \. 
When the above deterministic two-component system reaction-diffusion model is generalized to 
two dimensions, it leads to a Fisher genetic wave velocity v g = 2<\/D\a\, and a Fisher population 
wave velocities v = 2\J Da « v* = 2\/Da k . Assuming these waves are configured approximately 
at right angles as in Fig. [9J we obtain a y^-dependence for the bias m± = v g /v — y/a/a in the 
diffusion of the sector size, similar to Eq. (l47l) . In analogy with Eq. (1471) in the main text (with 
c— 1/2), we set a = sa, thus defining a selective advantage s under crowded conditions such that 
a a* . 

In contrast to Eq. (l47lh a linear dependence m± oc s can result if the competition between 
the mutant and wi l d typ e is weak, s < 1. In this limit, selection is weak compared to genetic 



drift iDoering et al.1 (120031 ). and the wave speed of a genetic Fisher wave is linear in s, 

v g = 2saDr g cid/ A , (weak selection) (BIO) 



where D is the usual spatial diffusivity of the organisms, c\d is the effective one-dimensional pop- 
ulation density at the frontier and r g is the generation time. The parameter A is the variance in 
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offspring number in a single generation, and is in all breeding models a number of order one. Note 
that which relation between m± and the selection coefficients (s or s) is realized strongly influences 
the genetic load, as predicted by Eq. f[58l) . 

Fisher population waves and Fisher genetic waves are approximately at right angles in the 
bacterial and yeast populations of Fig. [TJ Here, v « v* (so the population fronts advance into virgin 
territory at a common velocity) because these strains were chosen to be genetically neutral. There 
is no competition under crowded conditions, 5 = 0, and the Fisher genetic waves in this Figure are 
stalled out on average. Note that the model Eq. fjBID and its generalization to two dimensions do 
not apply to regions far behind the population front of the microbiological experiments of Fig. 1., 
because the used micro-organisms not only stop growing, but also stop diffusing once the population 
wave has passed. 



C. Simulations of genetic load at expanding frontiers 

In this section, we describe the simulations that were used to map out the genetic load at 
expanding frontiers as a function of the mutation rate, which is reported in Fig. [121 Our computer 
model is a derivative of the s o-called cont act process, which is believed to belong to the universality 



class of directed percolation lOdorl (120041 ) 



The simulation uses random sequential updates to evolve a one-dimensional lattice of binary 
variables cr^, i = 1 . . . TV, which can have the values WT (wild type) and MT (mutant). Initially, 
all sites are WT. At each time step, a pair of sites {c^, cr^ + i} is chosen such that periodic boundary 
conditions are respected. This pair of sites is then updated according to {c^,c^ + i} — > {c^, cr f i+1 } 
with certain transition rates w{a[ 1 c^ +1 |c^, cr^ + i). To describe the random occurrence of deleterious 
mutations, we implement 

w(MT,.\WT,.)=w(.,MT\.,WT)=2D x » d , (CI) 

where the states not shown are arbitrary. To account for the biased random walk of domain 
boundaries, 

w(WT, WT\WT,MT) = w(WT,WT\MT,WT) = (1 + m±)/2 (C2) 
w(MT, MT\WT, MT) = w(MT, MT\MT, WT) = (1 - m±)/2 . (C3) 

All other rates are zero. The expansion of domain boundaries is described by the first equation 
Eq. (IC2I) . and the shrinkage by the second Eq. (IC3I) . One lattice site in the simulation represents 
a length 2Dx of the linear frontier, and one time step corresponds to an equal change Ar = 2Dx 
in the frontier position, or effective time. The computer model then has the two non-dimensional 
parameters, 2Dx&d an d m±. 



We explore this model in the weak selection regime where the drift parameter m± <C 1 is 
small. This is the regime where the computer model generates domain boundaries that are well 
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approximated by slightly biased continuous time random walks, as assumed in the main part of 
the paper. We find that the error threshold, at which the wild type is lost, only depends on the 
parameter 7 = Dx^d/(^ m< j_) introduced in Eq. ([551) rather than on the two independent parameter 
of the computer model separately. This is documented by the data collapse in Fig. fT2l 
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